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ABSTRACT 

We present model independent determination of the CMB from any kind of fluid that has an 
equation of state taking four different values. The first region has w = 1/3, the second w = 1, 
the third w = —1 while the last one has —1 < w = cte < —2/3. This kind of dynamical w 
contains as a limit the cosmological constant and tracker models. 

We derive the model independent evolution of w^, for scalar fields, and we see that it remains 
most of the time in either of its three extremal values given by w$ = 1, — 1 ,u>t r . This ’’varying” 
w is the generic behavior of scalar fields, quintessence, and we determine the size of the different 
regions by solving the dynamical equations in a model independent way. 

The dynamical w models have a better fit to CMB data then the cosmological constant and the 
tracker models. We determine the effect of having the first two regions w = 1/3,1 and depending 
on the size of these periods they can be observed in the CMB. 

These models can be thought as arising after a late time phase transition where the scalar 
potential is produced. Before this time all the fields in this sector were massless and redshifted 
as radiation, giving the first period w = 1/3. 

In general, the CMB spectrum sets a lower limit to AN? and to the phase transition scale A c . 
For smaller AN? the CMB peaks are moved to the right of the spectrum and the hight increases 
considerably. 

Depending on the initial energy density we obtain a lower limit to the phase transition scale 
A c , when the scalar field appears and we have the transition from w = 1/3 to w = 1. For 
= 0.1 the CMB sets a lower limit to the phase transition scale A c > 0.2eV. For inverse 
power low potentials with <0.1 the constrain w a < —2/3 requires a power n < 1.8 and a 
phase transition A c > 4 MeV leaving a small energy scale window for models to work. 
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1 INTRODUCTION 


In recent time the cosmological observations on the cosmic microwave background radiation 
(”CMB”) [l| and the supernova project SNla pj] have lead to conclude that the universe is flat 
and it is expanding with an accelerating velocity. These conclusions show that the universe is 
now dominated by a energy density with negative pressure with = 0.7 ± 0.1 and w& < —2/3 
j|]. This energy is generically called the dark energy. Structure formation also favors a non¬ 
vanishing dark energy Q. 

It is not clear yet what this dark energy is. It could be a cosmological constant, quintessence 
(scalar field with gravitationally interaction) [9] or some other kind of exotic energy. 

The best way to determine what kind of energy is the dark energy is trough the equation of 
state parameter = 7^ — 1 = p/p, where p is the pressure and p the energy density of the fluid, 
and through its imprint on the CMB. The survey of redshifts of the different objects should in 
principle allow as to determine the value of w^o (o subscript referees to present day quantities) 
but only at small redshifts z. The result from the SN1A project pj] sets un upper limit to 
W(j )0 < —2/3 but does not distinguish a cosmological constant with (constant) w 0 = —1 and 
quintessence or any other form of matter with w 0 O < —2/3. It would be very interesting if in 
the future the SNla survey could constrain better the value of w ^ 0 . 

The CMB could give us information not only on the value of w^o but also on its form during all 
matter domination era. We will study models that have a changing w^ over time, well defined 
by the dynamics. We would like to see if we can distinguish from the CMB spectrum between 
a varying w^, a fixed w$ = cte and a true cosmological constant. Some general approaches can 
be found in Q. 

We will analyze the contribution to the CMB from a dark energy with a 7 ^ = w$ + 1 that takes 
four different values. It will have a w^ = 1/3 for energies above a certain scale A c , which we will 
call the phase transition scale. Starting at A c we will have a region with w$ = 1 and duration 
AW) where N is the logarithm of the scale factor a (N = Log[a\). Thirdly we will have = — 1 
for almost the same amount of time as in the previous period, AW — AW) and finally we will 
end up in a region with — 1 < w<p 0 = cte < —2/3 for a duration of AW- The cosmological 
evolution and the resulting CMB will have only four new parameters AW) AJV 2 , AW and w<p 0 . 
By varying these parameters we will cover a wide range of models. In particular we will cover 
all quintessence models. 

The analysis of the CMB with this kind of dark energy does not depend on its nature, it could 
be a scalar field (quintessence) or any other form of dark energy that gives the four sectors 
described above. However, we would like to point out that this pattern of different w$ is 
precisely what one expects from a quintessence scalar held and we will prove it. We will also 
show that it is generic, i.e. model independent. In the case of a quintessence scalar held the 
parameters AN\, AN 2 , AN^wso will be functions of A c , the phase transitions scale where the 
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scalar field is produced, the initial energy density of the scalar field, the minimum value 
of Umin = U(</>)/ 3H 2 and on the final value u^ 0 . For inverse power low potentials ”IPL” the 
number of parameters is reduced to three, A c , w^g. 


The evolution of scalar fields has been widely studied and some general approaches con be found 
in PH- The evolution of the scalar field (f> depends on the functional form of its potential 
V(4>) and a late time accelerating universe constrains the form of the potential [12]. Even though 
the evolution of the scalar field depends on the potential we will show that it is possible to obtain 
a model independent behavior of II 0 and w^. 


The contribution of scalar fields to the CMB has been studied but in most cases a constant w 0 
has been used. The fields with constant w 0 are called tracker fields Even though tracker 
fields are very interesting, specially because they do hardly depend on the initial conditions, 
they are not consistent with the observed (at least for inverse power potentials) . This work 
generalizes the tracker analysis since it contains the tracker model as a limiting case 


IPL tracker fields with constant w 0 = w^r = cte are not consistent with present day cosmological 
observations. Tracker fields require N > 5 and a small < —2/3 today requires n < 1. 
However, for n < 1 the scalar field has not reached its tracker value by present day. Of course, 
tracker fields are not the generic evolution of scalar fields. 


We will show that the generic behavior of 70 = w 0 + 1 for a quintessence scalar field with an 
arbitrary potential (with the restriction V > 0 and A i = —V'/V 3 > 1) has three critical points 
given by 70 = 2 , 0 and \ 2 Q ( p/3 (or w 0 = 1 , —1 and A 2 O 0 /3 — 1 ). 

The parameter 7 ^ will be most of the time in either of the three critical points. Independent 
of its initial value it will go rapidly to 7 ^ = 2 and remain there for a long period of time A N\. 
Afterwards it will sharply go to 7 ^ = 0 and stay there during almost the same amount of time 
as in the first stage AIV 2 — AN\. The amount of time it spends in these two regions depends 
only on the initial energy density, and on A c the phase transition scale. 

Finally, 7 ^ will evolve to its tracker value 7 ^ = \ 2 Q ( f ) /3 where it will remain. The amount of 
time before we reach present day, denoted by A N a , depends on the values of 0 ^ and 70 t r . 

Using this generic evolution of 7 ^ we can determine which models have the best fit to the acoustic 
CMB pattern by varying AN\, AN 2 , A N a and 70 O . The change of these four parameters covers 
all scalar field models. We hope to be able to infer form the results the phase transition scale 

A c . 

The work is organized as follows. In sectJH we give an overview of the models and we summarize 
the main theoretical and phenomenological results. In sect.|| we set the general dynamical 
equations for the quintessence field. In sect]4| we first analyze in a model independent way the 
evolution of 4> and then we do a model independent analysis of the dynamics of 7 ^ = + 1 . 

In both cases we give the model dependent parameters. In sect we show how long the regions 
with 70 = 2 ,O, 70 o l as t while in sect.||| we compare the CMB obtained in the presence of the 
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scalar field with the experimental data and with a true cosmological constant. Finally we give 
in sect.fj] our conclusions. 

If the reader is only interested on the setting of the model and its cosmological consequences 
she/he could skip sections || and |I| where we prove the model independence behavior of the 
70 = + 1 and go directly to sect.(|j). 

2 Overview 

The strategy is to analyze the spectra of CMB, using a modified version of CAMB, in a model 
independent way and see from its result if we can distinguish between different quintessence 
models, tracker, cosmological constant or other kinds of exotic energy densities. The results on 
the effect on CMB by the fluid with a generic behavior of 7 ^ = + 1, as seen from fig. will 

be valid independently of the nature of this fluid, i.e. scalar field or a exotic type of fluid. Notice 
that the w = 1/3,1,—1 ,wt r = cte model fits well with the numerical result of a IPL potential 
with n = 1 and = 0.05. 

In the case of a scalar held, we will assume that the scalar held appears at a scale A c with an 
energy density Q^(A C ). The late time appearance of the 0 held suggests that a phase transition 
takes place creating the scalar held. We are not concern with the precise mechanism of its 
appearance (see fn], |h|). However, energy conservation would suggest that the energy density 
of the <f> held after the phase transition would be given in terms of the energy density of the 
system before the phase transition and we will take them to be equal. It is natural to assume 
that all the energy density before the phase transition, in this sector, was in relativistic degrees 
of freedom. If the phase transition takes place after nucleosynthesis ”NS” then the primordial 
creation of nuclei puts un upper limit to the relativistic energy density to be less than 0 . 1 - 0.2 
of the critical energy density 0 .H- IfA c is larger than the NS scale then we do not need 
to worry about the NS bound since independent of its initial value, will drop rapidly and 
remain small for a long period of time (covering NS). 

In a chronological order, we would start with a universe filled with the SM particles and a Q 
sector (could be another gauge group) and with gravitational interaction between the two sectors 
only. In both sectors all fields start massless, i.e. they redshift as radiation. The evolution of 
the SM is the standard one and we have nothing new to say. However, the Q sector will have a 
phase transition at A c leading to the appearance of a scalar field (j> with a potential V(cf>), the 
quintessence held. Above A c the fields in this sector will behave as radiation. The evolution of <f> 
for energies below A c is that of a scalar held with given potential V. However, the precise form 
of V is unknown. In table [2] we show the different model independent regions that we consider. 
The model dependence lies only on the size of these different periods and on the value of 70 in 
the last region. 
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Sector 

Energy 

Duration 

7</> = + 1 

Radiation 

E> A c 

W> W 

4/3 

First 

E\ < A c 

N t < N < W 

2 

Second 

e 2 < w 

N\ < N < N 2 

0 

Third 

W < e 2 

N q < N < N 2 

A 2 fV3 


Table 1: We show the different regions, its duration and the value of 7 ^ in each region with 
N = Log[a\. With AW = N\ — N t and A W = W — W one obtains AW = dAN 2 , d=l,2 for 
matter or radiation dominance, and all model dependence is then given by AW: AW = N 0 — W 
and 70 in the third region (W is at present day). 


If a late time phase transition takes place, so that most of the time the universe has been 
dominated by matter, then AW — AW as seen from eq.(31). This will be the case for a 
transition scale A c smaller than the radiation-matter equality energy If A c E rm we 

have a large radiation domination epoch and then AW — 2AW- 


From a cosmological point of view we have only 4 free parameters AW: AW: AW and 7 for 
(the value of 7 ^ during the third period). With these parameters we cover all models. 

The cosmological parameters are given in terms of the field theoretical model dependent param¬ 
eters. 


3 Cosmological Evolution of o 


We will now determine the cosmological evolution of a scalar field 4> with arbitrary potential V(cj)) 
and with only gravitational interaction with all other fields. This field is called quintessence. 

The cosmological evolution of cj) with an arbitrary potential V((f > ) can be determined from a 
system of differential equations describing a spatially flat Friedmann-Robertson-Walker universe 
in the presence of a barotropic fluid energy density pb that can be either radiation or matter. 
The equations are 


H — --{pb+Pb + ft 2 ), 

p = -3 H(p + p), (1) 



where H is the Hubble parameter, cj) = dcfr/dt , p ( p ) is the total energy density (pressure). We 
use the change of variables x = and y = and equations 0 take the following form 
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Figure 1: We show the evolution of w^ and solid and dashed lines respectively for an IPL 
potential with n = 1 and ilpt = 0.05 in a matter background as a function of N = Log[a\. The 
dotted line represents the theoretical and we see that it makes a good fit to the numerical 
solution. Ni is given at the initial scale A c and N\, give the end of the regions with w<p = —1,1 
respectively while the solid vertical line at N a denotes present day. Notice that for N < Ni we 
are assuming that the energy density p^ redshifts as radiation and we are also assuming that 
radiation dominates for N < N^. 
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H 0]: 


x N = -3x ++ |x[2x 2 + 7 6 ( 1 -a: 2 -y 2 )] 
lIN = -]j^Xxy+ ^y[2x 2 + ^(1 - x 2 - y 2 )\ ( 2 ) 

H n = ~H[2x 2 + 7&(1 - x 2 - y 2 )] 

where N is the logarithm of the scale factor a, N = ln(a ); /at = df/dN for / = x,y,H ; 
7 b = 1 + rcfe and A (IV) = —V'/V with V’ = dV/dcj). In terms of x, y the energy density parameter 
is = x 2 + y 2 while the equation of state parameter is given by 7 ^ — 1 = = p^/p# = fryp' • 

It is clear that x 2 ,y 2 < 1. 

The Friedmann or constraint equation for a flat universe Q b + = 1 must supplement equations 
(| 2 j) which are valid for any scalar potential as long as the interaction between the scalar field and 
matter or radiation is gravitational only. This set of differential equations is non-linear and for 
most cases has no analytical solutions. A general analysis for arbitrary potentials is performed 
in 0 IH- All model dependence falls on two quantities: A(IV) and the constant parameter 
7 h = 1,4/3 for matter or radiation, respectively. We will be interested in studying scalar fields 
that lead to a late time accelerated universe, i.e. to quintessence, and in this case we will have 
a decreasing A(IV) jjl2| and a late time behavior A(IV) —► 0. For constant A(IV) (exponential 

potential) one can have an accelerating universe if A (IV) < y/6 but its dynamics would lead to 

an accelerating universe too rapidly, i.e. not at a late time as ours, unless we fine tune the initial 
conditions. 

It is also useful to have the evolution of = p^/ZH 2 = x 2 + y 2 and 7 ^ = l + w^ = 2cc 2 /(a ; 2 H-y 2 ), 
derived from eq. (I), @ 

(Q^)tv = 3(7 6 - - f^) (3) 

(7 <p)n = 37 ^( 7 ^ - 2) ^A^/|£ - 1 j (4) 


The evolution for the energy density, valid only for constant 7 ^, is the usual one 


P<j> 


P^ N ~ N ^- 


(5) 


and the evolution of f l<f, when it is much smaller than one and with constant 7 ^ is 


( 6 ) 


From now the subscript i stands for initial conditions, when the potential V appears, and the 
subscript o for present day values. 
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4 Model Independent Analysis 


4.1 Evolution of x,y and H 


We are interested in studying scalar potentials that lead to quintessence, i.e. a late time (i.e. 
present day) acceleration period of the universe. For this to happen one needs A = —m p iV'/V —> 
0 in the asymptotic limit (or to a constant less then one). An accelerating universe (slow role 
conditions) requires |A| < 1 and we want this period to be at a late time. We will consider 
potentials with V > 0 and since the <fi field evolves to its minimum V' < 0 and A > 0 where we 
are assuming, without loss of generality, models with cf> > 0. 

We will define the phase transition scale A c in terms of the potential by 

A c = m) V4 (7) 


where V) is the initial value of the potential and we will consider models that have an initial 
value 


A i — TTlpl 


V^i) 


> 1 . 


( 8 ) 


From dimensional analysis we expect A* = 0(m p i/A c ) 1. If we have a phase transition at a 

scale A c which leads to the appearance of the 4> field (e.g. composite field) then we would also 
expect (j)i ~ A c . We will be working with late time phase transition but A c could be as large as 
10 16 GeV and we will still have Aj 3> 1. 

An interesting general property of these models is the presence of a many e-folds scaling period 
in which A is practically constant and -C 1. 


A semi-analytic approach [20] is useful to study some properties of the differential equation 
system given by eqs. (!)• To do this we initially consider only the terms that are proportional to 
A, since A,; 1, then we follow the evolution of x, y and H so every period has a characteristic 

set of simplified differential equations. We see from eqs.(|2|) that the leading terms in x and y, 
for A 1, are xn = J §Ay 2 and yjv = — Xxy. Combining these equations we have 


x N x= VnU 


(9) 


with a constant circular solution 


n^ = x 2 + y 2 = x 2 (Ni) + y 2 (Ni ) = %(iV 4 ). (10) 

Since xjy is positive x will grow while y n is negative giving a decreasing y. This period ends at 
a scale N m i n with x 2 (N m i n ) ~ H^j(W) Since A i ^$> 1, the x and y derivatives are quite 

large and the amount of e-folds between the initial time with y* until y reaches its minimal value 
Vmin is very short. An easy estimate can be derived from 2 Mr/y = —cA S> 1, c = a/3/2x giving 
1 » Nmin — Ni = Log[y m in/yi\/c\i = 0(1/Xi), in the assumption cA* = cte. 
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The minimal value of y, given at N m i n , can be obtained from eq. dD with yN = 0. At his point 

we have _ _ 

H n _ ^ | 3~[76 + 1^(2 — 7 fc)] 1 

3 Hx 


A (Nr, 


( 11 ) 

2 < 1. We 


yA ^<jn l<pi 

where we have taken x 2 (N m i n ) ~ and Hn/H = — 3 / 2 ( 7 ;, + 11^(2 — 7 ;,)) 
see that A in eq. (0 is of order l/y/TI^i and we have Xi/X{N m in ) 1 . 

The value of y m j n depends on the functional form of V((j)), which sets the functional form of 
A = —V'/V. In general we have y m in = V {‘fimin) / 3 / HN n but without specifying V(cj)) it is not 
possible to determine y m in- 

For an inverse power law potential with V = A/ +n </> _?l = 2>y 2 H 2 one has 


Vmin — 


4+re /0 

A 2 x~ n / 2 
^ c rmin 


( 12 ) 


= Vi 




= Vi 


1 


where we have approximated ~ Hf = Vi/3yf = A^~ n (j) i n /3yf in eq. dD since N min -Ni < 
1 and we have taken from eq.([ll]) 4> m in = n/X m in — and 0, = n/Aj. If we assume that 

the initial value of 4>i = n/A* = nA c then eq.( 12 ) gives 


A, 


n/2 


Vmin — Vi 




(13) 


We see that ymin = 0{X i n - 2 ) ~ 0(A/^ 2 ) <C y, if is not too small. 

Shortly after y reaches its minimum value the scaling period begins. In this period we neglect 
the quadratic x,y <C 1 terms in eqs.(B) to find: 


VN 

y 


Hn 

' N 


(14) 


which leads to yH = H m i n ymin = cte. Notice that a constant Hy leads to a constant potential 
since V = 3 H 2 y 2 and therefore A and cj) will be constant during this scaling period, i.e. 


A (Nmin) - X(N 2 ) 


(15) 


where we have defined the scale N 2 as the end of the scaling period. Furthermore, still neglecting 
the quadratic terms on x and y in the third equation of system i) we get the expressions 


H = H min e~h»( N - N ^ 
y = ymine^ b{N ~ Nmin) ■ 


( 16 ) 


We can take in eqs. (HD Nmin — Ni and H min ~ Hi as discussed above, but y m in < Vi- 









( 17 ) 


In the same approximation (x, y <C 1) the evolution of x is given by 

x N = (-3 + ^7 b )x 

and we have a decreasing x since x n < 0 for all values of jb- The evolution is 

x(N) = x min (N min )e^ 3 + §^)PV-iW ( 18 ) 

The scaling period finishes when eq.(|b|) is no longer valid and the first term in y^ of eqs. (I) 
cannot be neglected. This will happen when Ax is of order one and x will be of the same order 
of y, i.e. 70 will be significant larger then zero (say 70 ~ 0 . 1 ). 

At the end of the scaling period we have l/x 2 ~ A (IV 2 ) = A (N m i n ) and 

fy(iv 2 ) = v 2 {n 2 ) + x 2 (n 2 ) ~ \(N min y 2 ~ % ( 19 ) 

as seen from eq. (0. The value of ^(A^) depends on the initial and can be much smaller 
than one. This happens in general for tracker fields since the growing of fij, from ^(A^) — 
Q<j>i "C 1 after the end of the scaling period to fl 0 o = 0.7 gives enough time for 70 to grow from 
70 ~ 0 to its tracker value 7 ^ t r = A 2 Ll ( f > /n 2 . On the other hand if \{N 2 ) is of the order one then 
x(N 2 ) ~ 1 and N 2 ) ~ = 0.7 and there is not enough time to allow 70 to grow to its 

tracker value and one has 0 < 70 O < 7 ^ r . 

When the scaling period is over, A and the field cj) start to evolve again to the minimum of V. 
70 grows and reaches its tracker value and may or may not remain constant for long period of 
time. At the end the late time behavior has A —> 0 and f~ y 2 —> 1 with 70 —>■ 0. 


4.1.1 Parameters 


There are only four independent parameters that fix the cosmological evolution of the models 
from its initial value to present day. These parameters are , A c , y m i n and the value of 70 O 
today. All other quantities can be derived from them. 


The amount of e-folds between the initial time W at A c and N±, the scale where w goes from 
w = 1 to w = —1, is set by the condition x ~ y and both are mush smaller than one. We use 
the evolution of x, from eqs.(16) and ( 111 ) to get 


N 1 -N i = i Log 


Vmin 


( 20 ) 


were we have assumed AT* ~ N m i n . Eq.(f20|) is independent of 76 . We can take x m i n = 
Vi — \/^m and for an IPL model we have y mvn ~ y*(A c /y 7 !^)”/ 2 and eq.(20) gives 


71 

ANl = Ah - A, = -Log 
fa 


y/Ompi 


( 21 ) 
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The amount of e-folds between the initial time Ni at A c and the end of the scaling period 
N 2 is given by eqs. ( 0 , (0 and ( 0 ) with y = y 2 ~ x 2 = 1/A (N min ) ~ ^0/“ and y min ~ 
Di{A c /y/n^~) n / 2 (with Ni ~ N min ) giving 


N 2 - Ni = —Log 
376 


y 2 

Vmin. 


= Log {> / °^ 1 


376 


A, 


( 22 ) 


and we have taken in the last eq.( 22 ) j/j ~ 

Finally, the amount of time between the end of the scaling regime A 7 and present day , which 
we denote be N 0 , with = 0.7 can be approximated from eq. (|l 6 |) by 


'Vo 

= ^Log 

^(j)0 

L2/2 J 

3 

Q(f)i 


AN 0 = N 0 -N 2 = —Log 
376 


where we have taken ~ y 2 since y 2 > x 2 and 75 = 1 today. 


(23) 


Summing eqs.(U) and (|23|) we have 


N 0 -N i = ^-Log[A- n ^J 2 - l) ^ 0 \ (24) 

which gives the total scale A Nt = N a — Ni between the initial time at A c and present day. 

We see that the size of the different regions can be determine by the four parameters A c , ymin 
and w^q. 


4.2 Evolution of 7 ^ = + 1 

We have seen the evolution of x, y, H in the preceding subsection and we would like now to show 
how W(j, evolves in a general framework. The tracker solution is just a special case (or the late 
time evolution) of the general behavior of the scalar field shown here. 

The evolution of the equation of state parameter, 7 ^ = 1 + as given by eq.(||]) has a generic 
behavior for all scalar fields independent of its potential. We see that ( 7 ^)jv = 0 has three 
solutions, 70 = 2 , 0 and X 2 Ll(j,/3 (or w= 1, —1 and X 2 Ll c f ) /3 — 1). 

The parameter 7 ^ will be most of the time in either of the three critical points. Independent 
of its initial value it will go quite rapidly to 7 ^ = 1 and remain there for a long period of time. 
The fast increase in 7 ^ is because Aj 1. Afterwards it will sharply go to 7 ^ = 0 and stay there 
during almost the same amount of time as in the first stage. Finally it will go to its last period 
given by the tracker value 70 tr = X 2 Ll ( f ) /3 where it will remain. 

The first stage ( 7 ^ = 2) represents a scalar field whose kinetic energy density dominates (E 
V), it is called the kinetic region, and the energy density redshifts as p$ ~ a -6 = e~ 6N . 
The second period ( 7 ^ = 0) is valid when the potential energy dominates (E^ -C V) and, 
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therefore, the field remains constant and the energy density redshifts as a cosmological constant 
with ~ o° ~ cte. The last critical value gives a 7 '^ tr = X 2 Q ( j ) /‘& and this is called the 
tracker value and it is not completely constant. The energy density redshifts as a tracker field 
p ^ ~ a - 3 TV>t 7 - = e ~ 37 V 70 tr _ yy e w j^ denote the beginning of the kinetic term by Ni and the end 
by N\ . The second period ( 7 ^ = 0) starts at N\ and hnishes at IV 2 . 

In table |] we summarize the behavior of p 11 ^, 7 ^ for the different cases. 

Let us dehne the quantity 

<25) 

We see from eq.([f]) that the sign of 7 ^ depends if A > 1 or A < 1. For A > 1 we have 7 </,jv > 0 
and the value 7 ^ m ax = 2 or wp = 1 , which is the maximum value for 7 ^, is a stable point, i.e. 
as long as A > 1 the parameter 7 ^ will grow towards its maximum value and will stay at this 
point. For A < 1 we have 7 ^ < 0 and the value 7 </> = + 1 = 0 will be a stable point also. 


First Period, 7 ^ = + 1 = 2 


In this first period one has 7 ^ = 2 and the redshift of <f> is much faster then radiation or matter 
and f Iff, will decrease. We will have at the end of the period IV = N±, 7^ ~ 2, A = 0(1) and 
h^(IVi) = n/(l + 77) with 


P<j>{N 1 ) = p<t>{Ni) 3(jv 1 -JVi)(2-7 !) ) 
Pb{N{) Pb (Ni) 


(26) 


see 


At the initial time since Aj > 1 we have A > X^jQ^i/6 1 since 7 ^ < 2. From eq.(Q) we 
that the derivative of 7 ^ 1 and 7 will rapidly go to its maximum value 2 . 

The kinetic period must stop at some point since A which is proportional to 0^ will decrease 
as well and it will eventually become less than one and the sign of 'ysN will become negative. 
The value of A at the beginning of the scaling regime (which is when y reaches its minimum 
value) is A = A min y / ' 0 ?i (IV min )/ 37 ^(IV m j n ) = [y b + fi^(2 - 7 b )]/2 (c.f. eq.Q) which is already 
smaller than one. However, even though A < 1 and ( 7 < 0, the period of 7 ^ ~ 2 remains 
valid for a long period of time since 7 '^(N min ) = 2x 2 min /(x 2 min + y 2 min ) ~ 2(1 - y 2 mi J x 2 min ) <C 1 
for = f ^$> y 2 nin . So we expect 7 ^ to be close to two for a long period of time and it will 

drop to one only until y ~ x. This can be quite large depending on the value of y m in , he. it 
depends on A c and the initial conditions yt since y m in = 0(yiX i n ^ 2 ) ~ 0(yi Ag 2 ) as can be seen 
from eq.([jj). The evolution of y in the scaling period is given by eq.(|l 6 |). 

How many e-folds IVi — Ni has this period depends on the initial conditions and on the phase 
transition scale A r . 
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Second Period, 7 ^ = + 1 = 0 


The second stage starts when x ~ y. We are still in the scaling regime with yH = cte and since 
we have <C 1 , as seen from eq. (© , A <C 1. The quantity 7 ^ has been decreasing and it will 
arrive at its minimum possible value 70 = 0 or wp = — 1 . As long as A < 1 the value of 7 ^ ~ 0 
will remain constant and the will be constant during all this time, this is the second part of 
the scaling regime. The transition time between 7 ^ = 1.9 and 7 ^ = 0.1 is quite short, about 
AN = 1, so we do not need to take into account the transition period. 

Notice that even though 7 ^ ~ 0 it is not completely zero since A 2 = A 2 fl ^/370 < 1 implies 
that 0 < < 37 ^/A 2 . Since in this period p$ redshifts much slower than radiation or matter, 

will start to increase and A will eventually become larger then one again. This is the end 
of period two. During all this period we have, 7 ^ ~ 0, 4> ~ cte , A = A min and the evolution of 
^</>(N 2 ) = r 2 /(l + r 2 ) is given by 

_ /O W' _ 3(W 2 -JVi) 7b 

2 “ »(i\y “»(%) e ' ,27) 


Since during this period the field 4> remains constant the value of A(JV 2 ) ~ A (Ni) ~ A (N min ) is 
also constant. 


The second period ends (as the scaling period) when eq.(14) is no longer valid 
term in the equation y n of eqs.(|2|) cannot be neglected. This happens for x(N 2 ) ~ 
discussion below eq.(|l 6 |)). 


and the first 
A(A hy 1 (c.f. 


Third Period, 7 ^ = + 1 = A 2 f^/3 


The third period starts when 7 ^ is not too small (i.e. x is comparable to y and 7 ^ = 0(1/10). 
During all this time we have A > 1 again. However, in this case it will 7 ^ will not arrive at the 
maximum value 7 ^ = 2 since A is not large and it will stabilize at 

7 (j>tr = A 2 Q^/3. (28) 


and we will have P^( N 0 ) = r^/{l + r$) with 


= P<j>i N o) = P^i. N 2 ) -3(JV 0 -JV2)(7«tr-7b) 
3 “ Pb(N 0 ) p b (N 2 ) 


(29) 


If 7 <j>tr < 7 b then > 1. While A 2 !!^ remains constant we have the constant tracker value for 
70 or W'f,. A constant 7 ^ is possible when 1. However, at late times the attractor value 

will be 70 tr —> 0 and —» 1 since is constrained to be smaller than one and A —> 0. But, 
even for 7 ^t r not constant the evolution of 7 ^tr in eq.(28) is valid and the value generalizes the 
tracker behavior. 


For an inverse power law potential V = V 0 cf> n we have A = n/(f> and 70 tr = n 2 ^/3 (/ 2 which is 
the valued obtained by Pi , [pcfl . 
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5 Duration of the Periods 


In order to know the relative size of the different periods we can use eqs.(|26l) and ([27|). Let us 
define A N± = N\ — Ni and A N 2 = N 2 — N\, they give the amount of e-folds during the 7 ^ = 2 
and 74 , = 0 periods, respectively. Combining both eqs.(|26|) and (f27|) we have 


L2 _ P<j>(N 2 )pb{Ni ) _ e -3AiV 1 (2-7,,)+3AiV276 

n Pb(N 2 )p r f > (N i ) 


(30) 


If the exponent in eq.(30) is zero than we will have r 2 
Solving for AN 2 in eq.(30) we obtain 


7 which implies that Q,p(N 2 ) = f 


AN 2 = ANi 


7 b 


-1 + —Log 


37 b 


[2 

Ti. 


(31) 


If we use the result of quintessence evolution at the beginning and end of the scaling period 
Qij>(Af 2 ) = ll 0 i(iVj) given in eq.([H]) we have r 2 = ' 7 . For matter, 75 = 1, and eq.(|3l|) gives 
A N 2 = ANi while for radiation, 7 b = 4/3, and A N 2 = ANi/2. 


The universe has been dominated by matter for a period of N a — N rm ~ 8 , where N a stands for 
present day value and N rm for the scale at radiation-matter equivalence. 


Including the third period we have from eqs.(26), (27) and (| 


n = P<p(N 0 )P b m = 3AAr 1 (2-7 6 (Ar 1 ))+3AAf 2 7 b (Ar 2 )-3AAr 0 (7^ r .-l) = V _A -3AAT o ( 70tr -l) ( n ? \ 

r% ptiNjp+m n 1 j 


where we have assumed that the third period is already at the matter dominated epoch, 7 &(IV 0 ) = 
1. It is clear from eq.([32|) that the size A N 0 and the value of 70 t r will set the initial energy 
density 7 , Q^A/) = 7/(1 + 7 ) assuming that the final stage of period three is today N 0 and 
n^N 0 ) = = 0.7 which gives r 3 = 7/3. If we take in eq.(32) the equality f l^Ni) = 

which implies 7 = r 2 , then we can express 7 = and = 7 /( 1 + 7 ). Of course, 

on the other hand, if we know then we can determine AN 0 r y ( j }tr . 


As a function of AN? = AN 2 + AN 2 + A N a we can estimate the magnitude of the phase 
transition scale A c . From A c = V)J^ = (3 Hfyf) 1 ^ and using the approximation that <S 1 
during almost all the time between present day and initial time (at A c ) we have 


H, = H n 


0 3"/bAN T /2 


(33) 


giving a scale 


A c = 3 1/4 Vyjf 0 e 37tA7VT/4 . 


(34) 


The scale A c increases with larger A Nt- From eqs.(|2^) and (34) we can derive the order of 
magnitude for A c in terms of n and H a giving A c = H 0 ^ <y ‘ i+n ' > ■ If we know A Nt then we can 
determine A c and the power n for IPL models. 
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How long do the periods last depends on the models and by varying the size of AN 2 , A N a and 
7 ^tr we cover all models. 

If A N 0 = 0 and AN 2 > A N rm = N a — N rm then the model would be undistinguishable from 
a true cosmological constant 7 ^ = w$ + 1 = 0 since during all the matter domination era the 
equation of state would be 7 ^ = 0. If we have A N a > AN rm then the model reduces to tracker 
models with a constant 7 during all the matter domination era. So, our model contains the 
tracker and cosmological constant as limiting cases. 

More interesting is to see if we can determine the nature and scale of the dark energy. For this 
to happen a late time phase transition must take place such that A c is at AA+ < AN rm . 


6 Analysis of CMB spectra 

We will now analyze the generic behavior of a fluid with equation of state divided in four different 
regions with w = 1/3,1, —1 ,wt r - We will vary the sizes of the regions and we will determine 
the effect of having regions with w = 1/3 or w = 1 in contrast to a cosmological constant or a 
tracker field (with —1 < wt r = cte < —2/3). This analysis is valid for all kinds of fluids with the 
specific equation of state and it is also the generic behavior of scalar fields. We will compare 
to the model wt r = —0.82 which was found to be a better fit to CMB than a true cosmological 
constant {J. 

6.1 Effect of Radiation Period, w = 1/3 

The first section we have w = 1/3 and the fluid (scalar field) redshifts as radiation. As long as 
the fluid has w = 1/3 its energy density will remain the same compared to radiation. If during 
nucleosynthesis the fluid has w = 1/3 then the BBNS bound requires the f l^NS) < 0.1 — 0.2 

0 H- 

In fig.([|) we show the different CMB for w = 1/3,0,—1 for A Nt = N a — N = 9, AIV 2 = 
N 2 — N\ = 4.5 and A N a = N a — N 2 = 0. We have chosen AA+ = 9 because it is the smallest 
value satisfying the condition AIVi = AIV 2 , AN a = 0 and giving the correct CMB spectrum. We 
have taken w = 1,-1 for N{ < N < N\ and N± < N < N 2 = N a , respectively. 

We see that the first and second peaks are suppressed for w = 1/3 compared to w = — 1 while 
the third peak is enhanced. The positions of the first two peaks is basically the same and the 
position of the third peak is moved from 868 to 864 (0.4%), for w = 1/3,—1 respectively. For 
smaller A/, i.e. more distant from present day, the effect is suppressed. It is not surprising since 
the N t would be further way from energy-matter equality and its effect on CMB would be less 
important. 
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Figure 3: We show the effect of the kinetic era with A/ < N < N\ by varying w = 1,0, —1 
and ANt = 2AN - 2 = 9. The vertical axes is l{l + l)ci/2ir([iK 2 ). 

In fig. (3) we show the CMB for different values of w = 1,0,—1 during N, < N < N\ and 
take w = 1/3 for N < Ni while w = — 1 for N± < N 2 = N 0 . The effect of having a kinetic 
period enhances the first three peaks and shifts the spectrum to higher modes, i.e. higher l. 
The curve for w = 0 is indistinguishable from the w = — 1 one. The position and hight of 
the peaks are pi = (227, 5275), p 2 = (559, 2605), p 3 = (868,2240) for w = 1 while for w = - 1 
we have p\ = (224,5138), p 2 = (545,2310), P 3 = (832,2165) giving a percentage difference 
pi = (1.3%, 2.6%), p 2 = (2.5%, 12.7%),P 3 = (4.3%, 3.4%). We see that the largest discrepancy 
is the altitude of the second peak. 

The difference in hight and positions may in principle distinguish between a cosmological con¬ 
stant and a scalar field, or any fluid with the specific equation of state behavior. 

The effect of having w ^ — 1 during < N < Aq is not significant for w = 0 and would be 
even less for a tracker fields with wt r = cte < —2/3 but it is observable for w = 1/3. 
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6.3 Equal length periods, AN\ = N\ — Ni = N 2 — N\ = AN 2 , A N 0 = 

n 0 -n 2 = 0 

We have studied the case with A N± = AN 2 , A N 0 = 0. In fig.(£|) we show the behavior for 
different values of A Nt = 2 AN 2 = 6, 8,9,12,16. 

There is a lower limit of A Nt that gives an acceptable CMB spectrum. The lower limit is 
A Nt > 9. For smaller A Nt the peaks move to the right of the spectrum and the hight increases 
giving a spectrum not consistent with the CMB data. 

For larger AN? > 9 the curves tend to the cosmological constant. It is not surprising since for 

l „ A at o A i\t :i- __ i_l i- i_ „ l_x: :xi„ ,1 j „ „„ „„ xi^ „x 



1 


Figure 4: We show the effect on the CMB by varying the AW = AN 2 = 3,4,4.5,6.8 with 
A N a = 0. The vertical axes is l(l + 1)ci/2tt(/j i K 2 ). 
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6.4 Scaling condition Q^Ni) = ^(A^) = 0.1 

Following the discussion in sec.([hl]), we now that a scalar field will end up its scaling period with 
a equal to its starting value, i.e. ^(IVj) = ^(A^) = 0.1. We have taken this value of fsince 
for N > Ni the energy density behaved as radiation and we have to impose the nucleosynthesis 
bound on relativistic degrees of freedom < 0.1 — 0.2. Imposing this condition we have 

determined the evolution of the CMB for three different values of wt r = —1,—0.82,—0.7. We 
have chosen to analyze the wt r = —0.82 because it was found to be the best fit tracker model 
by |J. We have w = 1/3 for N < Ni, w = 1 for Ni < N < N\, w = — 1 for N\ < N < N 2 



1 


Figure 5: We show the effect on the CMB by varying the ANt = N a — Ni = 5, 7.1, 9.5,12.5, 
with the constrain Q 0 j(A/) = = 0.1 and w^tr = — 1 ,AIV 2 = 1.03 and we include the 

cosmological constant w = —1, for comparison. The vertical axes is l(l + l)q/27r(/iA' 2 ). 

In fig.(|),(|) and (@) we show the curves for different values of A/ with the restriction that 
Q(j>(Ni) = ^(A^) = 0.1 and for = — 1, —0.82, —0.7, respectively. 
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Figure 6: We show the effect on the CMB by varying the A Nt = N q — Ni = 5.3, 7, 7.2, 7.4, 
with the constrain = ^(A^) = 0.1 and = —0.82, AIV 2 = 1.25. We also include the 

tracker with constant w == 0.82 The vertical axes is 1(1 + l)ci/2ir(fiK ' 2 ). 

















In the case of w = —1 we have AN 0 = 1.03 and that the smallest acceptable model has 
A Nt = 8.5, AN 2 = 3.6, see fig.(||). The best model has ANt = 8 . 88 , N\ — N a = 3.7 and peaks 
pi = (224,5133), p 2 = (549,2363),p 3 = (840,2178). 

For w = —0.82, fig.(||), the best fit for tracker models, we have AN 0 = N 0 — N 2 = 1.25 and the 
minimum acceptable distance is ANt = Ni — N a = 7.19, AN 2 = 3.9. Smaller values of ANt 
give a spectrum with peaks too large and second and third peaks moved to the right (high l 
modes). For large ANt the spectrum tends to the tracker spectrum w = —0.82. 

The best model has ANt = N\ — N 0 = 7.6, AIV 2 = 3.1 with peaks and position p\ = 
(223, 5016), P 2 = (551, 2450),P 3 = (851,2102). We have compared the x 2 of the models and 
the ANt = 7.6, = —0.82 model has a better fit than the tracker model with constant 

wtr = —0.82, which was found to be the best tracker fit Q. We see that having a dynamical w, 
is not only more reasonable from a theoretical point of view but it fits the data better. 

Finally, we consider w = —0.7 for N > A^- In this case we have A N 0 = N 0 — N2 = 1.47 
and the minimum acceptable model has ANt = 6 . 8 , AN 2 = 3.6, while the best model has 
AN t = 7.3, A N 2 = 3.8 with peaks p x = ( 222 ,4954), p 2 = (550, 2422), p 3 = (853,2035). 

We see that in all three cases w = —1,—0.82,—0.7, with condition Q^(A 7 *) = 11^(A^) = 0.1 we 
have a minimum acceptable value of ANt and for smaller ANt the peaks move to the right of 
the spectrum and the hight of the peaks increases considerably. This conclusion is generic and 
sets a lower limit to ANt, the distance to the phase transition scale A c , or equivalently it sets 
a lower limit to A c . 


The smallest ANt is set by the largest acceptable w (here we have taken it to be w<p 0 = —0.7) 
giving in our case a ANt = 6.8 for = 0.1. This result puts a constraint on how late the 
phase transition can take place. In terms of the energy A c = = [Iipj3i7 ? 2 ] l//4 we can set a 

lower value for the transition scale. Using eq.(34) with = 0.1 andAAU 1 = 6.82 we get 


A c = = 2 x 10~ 10 GeU = 0.2eU 


(35) 


i.e. for models with a phase transition below eq.(35) the CMB will not agree with the observa¬ 
tions. This result is independent of the type of potential. 


Furthermore, we now that for inverse power potential there is un upper limit to A c coming by 
requiring that w < —2/3. The limiting value assuming Qm < 0.1, for V = A 4 :+n (j)~ n , is n < 1.8 
giving A c = 4 x 10 _3 GeU ~ ^-2/(4+n) ^ 'p] iere f or . e5 £ or jp^ potentials the only acceptable models 
have phase transition scale 


4 x 10“ 3 GeU > A c > 2 x 10' 10 GeU. 


(36) 
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7 Conclusions 


We have analyzed the CMB spectra for a fluid with an equation of state that takes different 
values. The values are w = 1/3,1,-1 ,Wf r for N in the regions 7V ? ; — Npi, anc k, N\ — A^ — 
N\. N a — N 2 , respectively. The results are independent of the type of fluid we have. The 
cosmological constant and the tracker models are special cases of our general set up. 

We have shown that the evolution of a scalar field, for any potential that leads to an acceler¬ 
ating universe at late times, has exactly the kind of behavior described above. It starts at the 
condensation scale A c and enters a period with w = 1, then it undergoes a period with w = — 1 
and finally ends up in a region with — 1 < w < —2/3. We have shown that the energy density 
at the end of the scaling period (end of w = — 1 region) has the same energy ratio as in the 
beginning, i.e. = Q<f>(N 2 ). The time it spends on the last region depends on the value of 

Q<f>(N 2 ) and on w during this time. Before the phase transition scale A c we are assuming that 
all particles were at thermal equilibrium and massless in the quintessence sector. At the phase 
transition scale A c the particles acquire a mass and a non trivial potential. 

We have shown that models with w = 1/3,1, —1 ,wt r have a better fit to the data then tracker 
or cosmological constant. Furthermore, we have determined the effect of the first two periods 
w = 1/3 and w = 1 and even though the effect is small it is nonetheless observable. 

In general, the CMB spectrum sets a lower limit to ANt, which implies a lower limit to the phase 
transition scale A c . For smaller ANt the CMB peaks are moved to the right of the spectrum 
and the hight increases considerably. 

For any f1^, the CMB sets a lower limit to the phase transition scale. In the case of 54^ (A/) = 0.1 
the limit is A c = 0.2eC for any scalar potential. We do not take 0,^ much larger because we 
should comply with the NS bound on relativistic degrees of freedom < 0.1 — 0.2. If we take 
“C 0.1 then the constrain on the phase transition scale will be less stringent since the effect 
of the scalar field is only relevant recently <C 1 during all the time before present time). 

For inverse power law potentials we can also set un upper limit to A c and for < 0.1 it gives 
an inverse power n < 1.8 and A c < 4 x 10 ^^GeV. In this class of potentials only models with 
4 x \0~ z GeV > A c > 2 x 10 ~ 8 GeV would give the correct w^ 0 and CMB spectrum. 
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